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The  kinetics  for  the  reactions  of  NOO+  ions  with  neutral  molecules  having  ionization  potentials 
(IPs)  from  9.27  to  15.58  eV  was  measured  in  a  selected  ion  flow  tube  at  298  K.  The  NOO+  ions  are 
produced  from  the  reaction  of  N3+O2  and  have  been  reacted  with  the  following:  NO,  C6F6,  CS2, 
CF3I,  C3F6,  OCS,  C2H6,  Xe,  S02,  03,  N20,  C02,  Kr,  CO,  D2,  and  N2.  Numerous  types  of  reactions 
were  observed  with  the  various  neutral  reagents,  including  production  of  NO+  (which  may  involve 
loss  of  an  O  from  the  ion  or  addition  of  O  to  the  neutral  reactant,  although  the  two  channels  could 
not  be  distinguished  here),  charge  transfer,  isomerization  of  NOO+  to  ONO+,  and  hydride 
abstraction.  High  level  theoretical  calculations  of  the  structures  and  energetics  of  the  various 
isomers,  electronic  states,  and  transition  states  of  NOO  and  NOO+  were  performed  to  better 
understand  the  observed  reactivity.  All  neutral  species  with  an  IP  1 1.18  eV  were  observed  to  react 
with  NOO+  in  part  by  charge  transfer.  Detailed  calculations  showed  that  the  recommended  adiabatic 
and  vertical  IPs  of  NOO  are  10.4  and  11.7  eV,  respectively,  at  the  MRCISD(Q)/AVQZ  level  of 
theory.  The  observed  experimental  limit  for  charge  transfer  of  11.18  eV  agreed  well  with  the 
energetics  of  the  final  products  obtained  from  theory  if  dissociation  of  the  neutral  metastable  product 
occurred,  i.e.,  the  products  were  X++[0(3P)  +  NO(2II)],  where  [0(3P)+N0(2fl)]  formed  via 
dissociation  of  metastable  NOO.  Charge  exchange  with  neutral  reagent  X  would,  therefore,  be 
exothermic  if  IP(X)<[IPad(NOO)-A£(0+No).NOo]=  ~  1 1-1  eV,  where  IPad(NOO)  is  the  adiabatic 
IP.  The  potential  energy  surface  for  the  reaction  of  NOO+  with  C2H6  was  also  calculated,  indicating 
that  two  pathways  for  formation  of  HN02+C2H5  exist.  ©  2006  American  Institute  of  Physics. 
[DOI:  10.1063/1.2177257] 


I.  INTRODUCTION 

The  reaction  of  N3  with  02  proceeds  by  two  channels, 
N;  +  02->NO+  +  (N20)  (la) 

-^NOj  +  N2,  (lb) 

where  N20  is  shown  in  parentheses  because  several  different 
neutral  fragments  can  be  formed  exothermically.  Several 
studies  of  reaction  (1)  have  been  performed,  a  number  of 
which  have  indicated  that  the  N02  product  has  the  isomeric 
form  NOO+.1-7  Our  groups  have  recently  studied  this  reac¬ 
tion  over  an  extended  temperature  range  showing  that  chan¬ 
nel  ( 1  b)  diminished  substantially  with  increasing  temperature 
from  120  to  1400  K. 6  In  that  study,  it  was  conclusively 
shown  that  the  N02  isomer  formed  was  not  the  ONO+  form 
made  by  ionizing  ONO. 

"■Also  at  the  Institute  for  Scientific  Research,  Boston  College,  Chestnut  Hill, 
MA.  Electronic  mail:  anthony.midey@hanscom.af.mil 
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d,Electronic  mail:  albert.viggiano@hanscom.af.mil 

“'Present  address:  ITAMP  Harvard-Smithsonian  Center  for  Astrophysics, 
Cambridge,  MA. 
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While  reaction  (1)  has  been  studied  extensively,  little 
research  has  been  performed  to  elucidate  the  ion  chemistry 
of  NOO+.  Matusoka  et  a  if  found  in  a  time  resolved  atmo¬ 
spheric  pressure  ionization  (TRAPI)  instrument  that  the  N02 
ions  formed  in  reaction  (1)  were  lost  in  collisions  with  N2. 
The  experiments  showed  that  the  NOO+  reaction  with  N2 
formed  NO+  (assumed  to  be  from  thermal  dissociation  of 
NOO+)  and  that  the  reaction  had  an  activation  energy  asso¬ 
ciated  with  it.7  On  the  other  hand,  Hiraoka4  showed  that 
these  ions  do  not  react  in  a  high-pressure  mass  spectrometer 
(HPMS).  The  latter  results  were  consistent  with  ab  initio 
calculations  of  the  successive  binding  energies  of  N2  ligands 
in  N02(N2)„  clusters  containing  ONO+.b  In  addition,  the  cal¬ 
culations  further  indicated  that  NOO+  was  not  stable.  In  the 
TRAPI  experiments,  it  was  further  shown  that  the  N02  pro¬ 
duced  in  reaction  (1)  reacted  with  C2H6  to  form  C2Hj, 
whereas  ONO+  ions  did  not  react  with  C2H6,  again  indicat¬ 
ing  that  the  NOO+  isomer  is  formed  in  reaction  (l).7 

All  of  the  above  studies  were  made  in  complicated  mix¬ 
tures  of  gases  involving  multiple  steps  in  the  observed  reac¬ 
tions.  However,  in  our  recent  study  of  reaction  (I)  in  a  se¬ 
lected  ion  flow  tube  (SIFT),6  we  could  better  control  the 
chemistry  to  produce  NOO+  than  in  the  previous  experi¬ 
ments.  N3  ions  were  injected  into  a  pure  He  buffer  where  02 
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was  added  in  substantial  quantity  through  an  upstream  inlet 
to  produce  NO+  and  NO,,  neither  of  which  reacted  with  02. 
After  reaction  (1)  was  driven  to  completion,  Xe  was  added 
downstream  to  distinguish  between  the  two  N02  isomers. 
We  found  that  the  N02  formed  in  reaction  (1)  reacted  with 
Xe  and  that  ONO+  produced  from  ONO  did  not.  This  con¬ 
clusively  showed  that  N3  reacts  with  02  to  form  N02  in  the 
NOO+  form.6 

In  the  present  study,  the  ion  chemistry  of  NOO+  is  ex¬ 
amined  in  more  detail  in  a  SIFT  at  298  K,  where  reaction  (1) 
was  used  to  create  NOO+  in  the  flow  tube.  The  following 
molecules  were  used  as  neutral  reagents,  listed  in  order  of 
increasing  ionization  potential  (IP):  NO,  C6F6,  CS2,  CF3I, 
C3F6,  OCS,  C2H6,  Xe,  S02,  03,  N20,  C02,  Kr,  CO,  D2,  and 
N2.  Using  molecules  with  IP  values  ranging  from 
9.27  to  15.58  eV  brackets  the  IP  for  NOO  and  gives  insight 
into  the  reactivity  of  NOO+.  Numerous  types  of  reactions  are 
observed  with  the  various  neutral  reagents,  including  produc¬ 
tion  of  NO+  (which  may  involve  loss  of  an  O  or  addition  of 
O  to  the  reactant  neutral,  although  the  two  channels  could 
not  be  distinguished  here),  charge  transfer,  isomerization  of 
NOO+  to  ONO+,  and  hydride  abstraction.  Based  on  the  di¬ 
versity  of  the  observed  chemistry,  the  measurements  have 
been  complemented  by  high  level  theoretical  calculations  of 
the  structures  and  energetics  of  the  various  isomers,  elec¬ 
tronic  states,  and  transition  states  of  NOO  and  NOO+.  In 
addition,  the  potential  energy  surface  for  the  reaction  of 
NOO+  with  C2H6  has  also  been  calculated. 

II.  EXPERIMENTAL  METHOD 

The  measurements  were  made  in  a  selected  ion  flow  tube 
(SIFT)  that  has  been  thoroughly  described  previously.9,10 
Briefly,  N3  ions  were  generated  from  N2  in  an  effusive 
source  in  a  high-pressure  source  region.  Two  mechanisms  for 
N3  production  were  operative:  three-body  association  of  N* 
with  N2  and  Penning  ionization  of  N2,  the  latter  of  which  is 
shown  in  Eq.  (2), 

N2  +  N2-^NJ  +  N.  (2) 

The  N3  ions  were  mass  selected  and  injected  into  the  flow 
tube  through  a  Venturi  inlet,  resulting  in  a  small  amount  of 
N+  from  dissociation  upon  injection  (<2%).  Helium  buffer 
gas  (AGA,  99.995%)  carried  the  ions  downstream,  where  02 
(Mass.  Oxygen,  99.999%)  was  added  82  cm  from  the  sam¬ 
pling  orifice.  The  02  converted  the  N3  ions  into  NO+  and 
NOO+  through  reaction  (I).  Sufficient  02  was  added  such 
that  the  reaction  was  driven  to  completion  before  either  of 
the  two  reactant  inlets,  which  were  located  at  36  and  52  cm 
from  the  sampling  orifice.  The  reactant  gas  was  added 
through  one  of  these  inlets  and  the  kinetics  was  mostly  stud¬ 
ied  in  the  normal  fashion,  since  neither  NO+  nor  02  gener¬ 
ally  reacted  with  the  ionic  products  or  neutral  reagents  stud¬ 
ied.  The  average  relative  error  in  the  rate  constants  is  ±15% 
and  the  absolute  error  is  ±25%. 10  All  of  the  neutral  reagents 
were  used  as  obtained  from  the  manufacturer. 

A  few  complications  did  arise.  In  the  reactions  of  NOO+ 
with  S02  and  CS2  some  ONO+  was  produced  through 
isomerization.  This  ion  obviously  has  the  same  mass  as 


FIG.  1.  Plot  of  total  ion  signal  at  m/z=46  amu  vs  concentration  of  S02  for 
the  reaction  of  N00++S02  measured  in  a  selected  ion  flow  tube  (SIFT)  at 
298  K. 

NOO+  (m/z= 46),  resulting  in  curvature  in  the  rate  constant 
plots.  An  example  for  the  S02  reaction  is  shown  in  Fig.  1. 
Since  ONO+  did  not  react  with  S02,  the  46  amu  signal  lev¬ 
eled  off  at  higher  S02  flows  when  all  of  the  NOO+  was 
reacted  away.  Thus,  the  ratio  of  the  N02  signal  in  the  flat 
portion  of  the  curve  to  that  portion  with  no  S02  added  gave 
the  branching  fraction  for  ONO+  production.  To  obtain  the 
rate  constant  k,  the  N02  signal  in  Fig.  1  observed  after  reac¬ 
tion  time  t  was  fitted  to  Eq.  (3), 

[46  amu]  =  [46  amu]0{exp(-i[so2],)(l  -  %  ONOVlOO) 

+  (%ON07100)},  (3) 

which  also  yields  the  branching  into  the  two  product  chan¬ 
nels.  Another  complication  arose  in  the  Xe  reaction,  where  a 
small  amount  of  the  observed  Xe+  signal  was  attributed  to 
the  slow  reaction  of  02  (formed  from  N+  present  in  the  flow 
tube  reacting  with  02)  with  Xe.  However,  the  contribution  of 
the  02  reaction  to  the  Xe+  signal  can  be  easily  subtracted. 
Lastly,  the  NO+  formed  in  a  number  of  the  NOO+  reactions 
had  a  background  contribution  due  to  reaction  (1).  Conse¬ 
quently,  the  NO+  formed  in  the  reactions  with  NOO+  is  sim¬ 
ply  the  difference  between  the  amount  of  NO+  observed  with 
and  without  the  reactant  neutral  added. 

III.  CALCULATION  METHOD 

As  will  be  discussed  later,  the  ionization  potential  (IP)  of 
NOO  is  very  crucial  to  understanding  the  charge  transfer 
mechanism  and  one  of  the  main  goals  of  the  present  study  is 
to  determine  its  value.  Therefore,  we  calculated  the  IP  at 
various  levels  of  theory.  More  specifically,  the  geometries  of 
ground  state  NOO(2A')  and  NOO+(12+)  were  optimized  at 
the  CASPT2  and  ROHF-UCCSD(T)  levels  of  theory  along 
with  Dunning’s  correlation  consistent  polarized  valence 
quadruple-zeta  basis  set  augmented  with  diffuse  functions 
(aug-cc-pVQZ),  denoted  as  AVQZ.11  Here,  CASPT2  refers  to 
internally  contracted  second-order  multireference  perturba¬ 
tion  theory.1'  In  the  CASPT2  calculation,  all  the  orbitals  ex¬ 
cept  the  Is  were  involved  in  the  correlation  treatment  and  the 
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TABLE  I.  Kinetics  of  the  reactions  of  NOO+  with  neutral  molecule  X  at  298  K.  The  Su-Chesnavich  collision  rate  constant  k,.  (Refs.  28  and  29)  is  given  in 
brackets. 


X 

IP  (eV) 

k  [kcY 

(cm3  s  ') 

Products 

Branching 

ratio 

A  H™  K 

rxn 

(kJ  mol-1) 

X 

IP  (eV) 

k  [kJ 
(cm3  s  !) 

Products 

Branching 

ratio 

Atf™  K 

(kJ  mol  ]) 

NO 

8.9  [7.6] 

NO+  +  NOO 

>0.98 

-1 18.7b 

C2H6 

10  [12] 

C2H]+HN02 

1.00 

-499 

9.27 

Ot+N20 

<0.02 

-251 

11.50 

c6f6 

8.2  [12] 

QFJ+NOO 

0.91 

-48 

Xe 

5.7  [8.1] 

NO++XeO 

0.97 

-350c 

9.90 

no++c6f6+o 

0.09 

-211 

12.13 

XeO++NO 

0.03 

! 

O 

00 

CS2 

15  [13] 

no++cs2+o 

0.17 

-174 

S02 

10  [18] 

NO++(S03) 

0.55 

-174 

10.07 

SO++CS+NO 

0.14 

-155 

12.35 

ono++so2 

0.45 

-434 

ONS++(OCS) 

0.03 

-686d 

S02  +  (CN+S) 

0.11 

79c 

03 

5.1  [10] 

N0++202 

0.98 

-566 

CSt+NOO 

0.41 

-42 

12.53 

Of+NOj 

0.02 

-314 

ono++cs2 

0.14 

-434 

N20 

0.15  [9.0] 

no++n2+o2 

1.00 

-505 

CF3I 

11  [14] 

N0+h-CF3I+0 

0.55 

-174 

12.89 

no++n2o+o 

-174 

10.28 

CF]+ONO  +  I 

0.18 

-280 

OP+CF3  +  NO 

0.01 

— 198f 

C02 

<0.002 

No  reaction 

fi++cf2+ono 

0.12 

-38 

13.78 

CF,I++NOO 

0.14 

-27 

Kr 

<0.03 

No  reaction 

c3f6 

9.8  [14] 

no++c3f6+o 

0.73 

-174 

14.00 

10.60 

c2f4n++co2+f2 

0.10 

55s 

C3FJ+(FN02) 

0.04 

-109h 

CO 

7.9  [8.0] 

NO++CO+0 

0.95 

-174 

C3FJ+NOO 

0.13 

10' 

14.01 

OJ+NCO 

0.05 

-1 

OCS 

8.9  [13] 

OCS++NOO 

0.18 

61 

d2 

<0.02 

No  reaction 

11.18 

SO^+CO+NO 

0.23 

-289 

15.47 

NO++CO+SO 

0.59 

-391 

n2 

<0.008 

No  reaction 

15.58 


“All  rate  constants  shown  are  X  10-l°. 

bCharge  transfer  reaction  enthalpies  calculated  assuming  neutral  NOO  for  all  cases. 
cTaken  from  Ref.  6. 

dAWf(298  K)  ONS+  from  G3(B3LYP)  calculation. 
cA/ra298  K)  SO]  from  AH1}  S02  and  IP  of  S02. 
fA//y(298  K)  OI+  from  Hassanzadeh  and  Irikura  (Ref.  31). 
sA7/?(298  K)  C2F4N+  from  G3(B3LYP)  calculation. 
bA//7(298  K)  C3F5  upper  limit  from  Jarvis  et  at.  (Ref.  32). 

*A/^(298  K)  C3Fj  from  A fl®  C3F6  and  IP  of  C3F6. 


reference  wave  function  was  obtained  from  the  correspond-, 
ing  full  valence  state  specific  complete  active  space  self- 
consistent  fieldli'14  (CASSCF)  calculation,  which  included 
17  and  16  electrons  distributed  over  12  orbitals  for  the  neu¬ 
tral  and  ionic  species,  respectively.  ROHF-UCCSD(T)  refers 
to  the  spin  unrestricted  open  shell  coupled-cluster  singles 
and  doubles1'’  augmented  by  triple  excitation 
noniteratively.16  In  this  scheme,  the  standard  RHF-CCSD(T) 
is  used  for  the  closed  shell  systems.  For  the  open  shell  sys¬ 
tem,  the  ROHF  calculation  was  performed  first,  from  which 
unrestricted  molecular  orbitals  (MOs)  were  created,  followed 
by  UCCSD  calculations.  To  further  refine  the  energetics,  a 
single-point  multireference  configuration  interaction 
MRCISD(Q)  calculation  was  carried  out  on  top  of  the 
CASPT2  optimized  geometry  using  the  same  AVQZ  basis 
set.  MRCISD(<2)  refers  to  the  internally  contracted  multiref¬ 
erence  configuration  interaction  including  both  single  and 
double  excitations17,18  (MRCISD)  plus  the  multireference 


version  of  the  Davidson  correction  (Q).19  Similar  to  the 
CASPT2  calculation,  only  the  Is  orbitals  were  excluded  in 
the  correlation  treatment  and  the  same  full  valence  state  spe¬ 
cific  CASSCF  wave  functions  also  served  as  the  references. 
All  the  calculations  were  performed  using  the  MOLPRO  2002.6 
quantum  chemistry  program  package." 

IV.  RESULTS  AND  DISCUSSION 
A.  Experimental  results 

Rate  constants  and  product  branching  ratios  for  the 
NOO+  reactions  studied  at  298  K  are  listed  in  Table  I.  The 
neutral  reagents  are  listed  in  order  of  increasing  ionization 
potential  (IP).  The  heats  of  reaction  at  298  K  shown  in  Table 
I  have  been  calculated  using  the  standard  heats  of  formation 
available  in  the  NIST  Chemistry  Webbook,21  except  where 
noted.  For  the  less  common  species  C2F4N+  and  ONS+,  no 
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TABLE  II.  Optimized  NOO  and  NOO+  geometries  (bond  lengths  R  in  A  and  bond  angles  A  in  degrees)  at 
different  levels  of  theory.  In  all  the  multireference  methods,  full  valence  CASSCF  space  including  16  e~  (lor 
NOO+)  or  17  e~  (for  NOO)  in  12  orbitals  was  used  in  the  correlation  calculation. 


NOO(2A') 

NOO+('S+) 

Rft-O 

Ro-o 

^N-O-0 

Rn-o 

Ro-o 

•4n-o-o 

CASSCF/AVTZ 

1.2230 

1.3709 

120.3 

1.1387 

1.2321 

180.0 

CASPT2/AVQZ 

1.2315 

1.3482 

122.6 

1.1501 

1.2289 

180.0 

RHF-UCCSD(T)/AVQZ 

1.2128 

1.3451 

122.2 

1.1379 

1.2288 

180.0 

MRCISD((7)/cc-pVTZa 

1.1459 

1.2334 

180.0 

“Results  are  taken  from  Ref.  6. 


literature  values  for  the  heats  of  formation  were  found. 
Therefore,  the  standard  heats  of  formation  at  298  K  for  these 
species  were  calculated  using  Gaussian  03  for  optimized 
geometries.22  Both  ONS  and  ONS+  were  observed  experi¬ 
mentally  by  Andrews  et  al.  using  infrared  (IR)  spectroscopy 
and  they  were  shown  to  have  bent  and  linear  structures,  re¬ 
spectively,  as  verified  by  density  functional  theory  (DFT) 
calculations  at  the  B3LYP  level.22  Consequently,  G3(B3LYP) 
theory24  was  used  to  calculate  the  heat  of  formation  of 
ONS+,  yielding  geometries  for  ONS  and  ONS+  in  good 
agreement  with  the  previous  calculations.2’  Analogous 
G3(B3LYP)  calculations  were  done  for  C2F4N  and  C2F4N+, 
the  latter  of  which  had  been  observed  in  trace  amounts  in  the 
mass  spectra  from  the  analysis  of  perfluorocarbon 
compounds.25"27  The  minimum  energy  isomer  of  C2F4N+ 
was  located  by  first  performing  optimizations  of  1 1  possible 
isomers  of  the  neutral  and  the  cation  at  the  B3LYP/6-311 
+  +G(3c//,2p)IIB3LYP/6-31G(r/)  level.  The  wave  function 
stability  of  all  species  needed  to  calculate  the  thermochem¬ 
istry  was  verified  as  well.  The  heats  of  formation  at  298  K 
obtained  from  these  calculations  are  976  ±5  kJ  mol-1  for 
ONS+  and  185±5  kJmor1  for  C2F4N+. 

Most  of  the  NOO+  reactions  occurred  with  a  rate  con¬ 
stant  k  at  or  near  the  Su-Chesnavich  collision  rate  constant 
limit28,29  kc,  i.e.,  k/kc> 0.5  or  50%  efficient.  C02,  Kr,  D2, 
and  N2  were  found  to  be  unreactive.  N20  was  the  only  spe¬ 
cies  that  reacted  slowly,  with  a  rate  constant  around  1.7%  of 
the  collisional  value,  kc.  For  the  species  with  the  lowest  ion¬ 
ization  potentials  (NO  and  C6F6),  the  reaction  proceeded  pre¬ 
dominantly  via  nondissociative  charge  transfer  (>91%  of 
products).  For  species  with  ionization  potentials  between 
10.07  and  11.18  eV,  charge  transfer  was  still  important 
( 1 8%— 41  %  of  the  products).  For  species  with  IP  >  11.5  eV, 
no  charge  transfer  occurred.  This  observation  implies  that  the 
IP  of  NOO  is  <11.5  eV.  The  calculations  discussed  below 
address  this  issue  in  more  detail. 

Production  of  NO*  ions  was  found  in  the  reactions  of 
NOO+  with  NO,  C6F6,  CS2,  CF3I,  C3F6,  OCS,  Xe,  03,  N20, 
and  CO.  It  was  the  predominant  reaction  channel  for  all  spe¬ 
cies  with  IP  2=  10.28  eV,  except  for  C2H6,  where  no  NO+  was 
formed.  Dissociation  of  NOO*  into  NO*  plus  O  is  exother¬ 
mic,  but  a  substantial  barrier  to  unimolecular  dissociation  of 
NOO*  exists.6  However,  in  the  course  of  an  ion-molecule 
reaction,  the  barrier  might  disappear  and  the  NO+  product 
might  be  the  result  of  either  O  atom  transfer  or  simple  dis¬ 


sociation  induced  by  the  proximity  of  the  reactant  neutral. 
Our  experiment  could  not  differentiate  between  the  two  pos¬ 
sibilities. 

Several  of  the  reactions  were  quite  complicated,  result¬ 
ing  in  three  to  six  different  product  ions.  Several  unusual 
product  species  were  observed,  including  ONS+  (with  CS2), 
FI+  and  IO+  (with  CF3I),  C2F4N+  (with  C3F6),  and  XeO* 
(with  Xe).  The  reactions  of  CS2  and  S02  caused  isomeriza¬ 
tion  of  NOO+  to  ONO+  in  14%  and  45%  of  the  reactions, 
respectively. 

B.  Ionization  potential  calculation 

The  optimized  equilibrium  geometries  for  NOO  and 
NOO*  are  given  in  Table  H  and  the  calculated  vertical  and 
adiabatic  IPs  are  collected  in  Table  III.  In  general,  the  geom¬ 
etries  from  all  of  the  methods  used  here  agreed  well  with 
each  other.  The  most  significant  difference  was  the  N-0 
bond  length.  The  multireference  based  correlated  method 
yielded  a  longer  bond  length  (—0.012  A).  The  structure  of 
NOO*  had  been  investigated  in  our  previous  study  at  the 
MRCISD((2)  level  of  theory  with  the  valence  triple-zeta  ba¬ 
sis  set  (cc-pVTZ).6  The  optimized  N-0  distance  was 
1.146  A,  which  lies  in  between  the  results  of  the  CASPT2 
and  ROHF-UCCSD(T)  calculations. 

In  contrast,  the  calculated  IP  values  showed  relatively 
large  differences  amongst  the  different  theoretical  methods, 
especially  for  the  adiabatic  IP.  The  general  trend  was  that  the 
MRCISD(<2)  calculations  predicted  smaller  values  than  the 
other  two  methods.  The  difference  between  the  MRCISD((2) 
and  CASPT2  results  was  only  —0.1  eV.  However,  the  most 

TABLE  III.  Adiabatic  and  vertical  ionization  potentials  (IPs)  of  NOO(vi') 
calculated  at  different  levels  of  theory.  In  all  the  multireference  methods, 
full  valence  CASSCF  space  including  16  e~  (for  NOO+)  or  17  e~  (for  NOO) 
in  12  orbitals  was  used  in  the  correlation  calculation. 


Vertical  IP  (eV) 

Adiabatic  IP  (eV) 

w/ZPE’ 

w/o  ZPE 

w/ZPE  w/o  ZPE 

MRCISD/AVQZ//CASPT2/AV  QZ 

11.69 

11.75 

10.44 

10.35 

CASPT2/AVQZ 

11.75 

11.81 

10.53 

10.44 

ROHF-UCCSD(T)/AV  QZ 

11.97 

12.05 

11.46 

11.37 

“ZPE  was  only  calculated  at  CASSCF/aug-cc-pVTZ  level  of  theory  based  on 
the  same  CASSCF  method  optimized  geometry.  As  for  the  vertical  IP,  pro¬ 
jected  frequency  calculation  was  carried  out  also  at  the  CASSCF  level  to 
obtain  the  ZPE  value. 
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significant  discrepancy  was  for  the  adiabatic  IP  obtained 
from  the  ROHF-UCCSD(T)  calculation,  which  was  —1.0  eV 
higher  than  the  MRCISD(g)  results.  In  our  previous  inves¬ 
tigation  of  NOO*,  consistent  results  amongst  the  coupled- 
cluster,  CASPT2,  and  MRCISD(<2)  methods  were  obtained 
for  the  0-0  bond  cleavage  process.6  Therefore,  the  differ¬ 
ence  was  likely  from  the  description  for  the  neutral  open 
shell  NOO  species.  The  7)  value  '0  of  0.05  in  the  UCCSD(T) 
calculation  of  NOO  indicates  the  extensively  multireference 
nature  of  this  species,  which  most  likely  caused  the  adiabatic 
IP  to  be  overestimated.  Therefore,  the  suggested  values  for 
adiabatic  and  vertical  IPs  of  NOO  are  10.44  (10.35)  and 
11.69(11.75)  eV,  respectively,  from  the  MRCISD(<2)  calcu¬ 
lations,  where  the  values  in  parentheses  do  not  include  the 
zero-point  energy  (ZPE)  correction.  At  this  level  of  calcula¬ 
tion,  the  errors  are  expected  to  be  around  0. 1  eV. 

These  results  provide  a  qualitative  interpretation  for  the 
IP  dependence  observed  in  the  experimental  charge  transfer 
branching  ratios.  Charge  transfer  is  essentially  nonadiabatic. 
Thus,  the  multireference  methods  are  best  suited  to  study  this 
process.  However,  in  the  present  system,  it  is  too  computa¬ 
tionally  expensive  to  perform  comprehensive  investigations 
at  this  level  of  theory.  Consequently,  a  potential  energy  sur¬ 
face  (PES)  scan  was  performed  to  qualitatively  explain  the 
trends  observed  in  the  experiments. 

As  seen  in  Table  I,  the  fraction  of  reaction  occurring  by 
charge  transfer  decreased  with  increasing  IP  and  shut  off  at 
IPs  over  11.18  eV.  Using  the  reaction  of  NOO+  with  NO  as 
an  example,  we  scanned  the  PESs  involved  at  the  same 
MRCISD(<2)  level  of  theory.  The  scan  was  carried  out  with 
the  NOO  and  NO  species  at  infinite  separation.  Two  states 
corresponding  to  the  NOO+NO*  and  NOO*+NO  asymp¬ 
totes,  respectively,  were  scanned  along  the  NOO  bending 
coordinate  with  all  other  degrees  of  freedom  fixed  at  the 
equilibrium  geometries  of  NOO*('2*)  and  NO(2ll).  The  re¬ 
sults  are  shown  in  Fig.  2.  Clearly,  the  two  asymptotes  cross 
in  the  vicinity  of  the  linear  NOO*('2*)  minimum  region. 
Upon  bending,  these  two  states  will  cross  at  the  minimum  of 
the  seam  of  crossing  (MSX)  marked  in  Fig.  2.  The  energy  of 
this  MSX  is  about  —0.15  eV  relative  to  the  NOO+(‘S+) 
+NO(2ri)  asymptote,  which  is  below  the  ZPE  level 
(0.23  eV)  of  NOO*('2*).  During  the  approach  of  the  NO 
species,  NOO+  and  NO  might  form  a  weak  van  der  Waals 
complex  with  little  change  in  the  geometries  of  the  indi¬ 
vidual  species.  On  the  contrary,  this  region  corresponds  to  a 
high-energy  regime  of  the  PES  originating  from  NOO 
+NO+  because  of  the  nonlinear  minimum  of  the  NOO  spe¬ 
cies.  Thus,  this  crossing  point  will  be  sampled  every  time  the 
NOO+  molecule  vibrates.  At  infinite  separation,  however,  the 
nonadiabatic  coupling  element  between  these  two  manifolds 
is  zero.  As  the  two  species  approach,  the  degeneracy  of  these 
states  will  be  lifted  and  a  nonvanishing  coupling  element 
between  states  with  the  same  spin  symmetry  will  grow,  re¬ 
sulting  in  charge  exchange.  Besides  the  crossing  on  the  bend¬ 
ing  coordinate,  crossing  on  other  coordinates  (N-0  and  0-0 
stretchings)  should  provide  little  contribution  because  of  the 
higher  energy  requirements  as  shown  in  the  previous 
MRCISD(e)  study  of  NOO*.6 

It  should  be  noted  that  the  IP  of  NO  is  9.27  eV,  which  is 


NOO  Bending  Angle  (Degree) 


FIG.  2.  Potential  energy  vs  NOO  bending  angle  for  the  NOO+  +  NO  reaction 
system  at  the  MRCISD((2)  level  of  theory  (see  text  for  details). 

lower  than  the  adiabatic  IP  of  NOO.  Thus,  the  crossing  be¬ 
tween  these  two  PESs  must  be  higher  in  energy  than  the 
asymptote  of  NOO*+NO,  as  schematically  shown  in  Fig.  3. 
The  charge  transfer  minimum  seam  of  crossing  (CT-MSX)  is 
located  in  the  vicinity  of  the  NOO++NO  potential  well,  rep¬ 
resented  more  generally  as  NOO++X  in  Fig.  3.  As  the  IP 
value  for  neutral  molecule  X  increases,  the  second  curve  rep¬ 
resenting  the  PES  for  NOO+X*  (curve  2)  will  begin  to  shift 
higher  in  energy,  as  illustrated  by  the  dashed  line.  As  a  con¬ 
sequence,  a  charge  transfer  barrier  will  arise  and  charge 
transfer  should  cease  to  be  important  at  neutral  IP  values 
above  the  adiabatic  IP  of  NOO.  However,  this  mechanism 
seems  to  contradict  the  experimental  results  that  showed 
charge  transfer  for  several  species  with  adiabatic  IPs  above 
that  of  NOO. 


FIG.  3.  Schematic  potential  energy  surfaces  for  the  reaction  of  NOO++X, 
labeled  as  1.  Solid  line  curve  labeled  as  2  is  for  NOO+NO+  and  the  dashed 
curve  labeled  as  2'  is  for  X  with  a  much  higher  ionization  potential  (IP). 
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FIG.  4.  ZPE  corrected  potential  energy  surface  for  the  dissociation  of  NOO 
at  MRC1SD(Q)/AVQZIICASPT2/AVQZ  level  (energetics  in  parentheses 
are  not  ZPE  corrected).  Numbers  on  the  structures  refer  to  the  bond  length 
(A)  and  bond  angles  (degree).  The  imaginary  frequencies  at  TSs  are  given  in 
the  parentheses. 

This  observation  subsequently  implies  that  another 
mechanism  must  be  occurring,  at  least  for  molecules  with 
IPs  above,  or  even  near,  the  adiabatic  IP  of  NOO.  The 
mechanism  outlined  in  Fig.  3  assumed  that  NOO  was 
formed.  However,  NOO  is  metastable,  as  shown  in  the  po¬ 
tential  energy  diagram  in  Fig.  4.  There  are  two  dissociation 
processes  associated  with  the  ground  state  NOO,  the  N-0 
and  0-0  bonding  fissions.  Two  corresponding  transition 
states,  TS1  and  TS2,  were  located  at  the  CASPT2/AVQZ 
level  of  theory  as  shown  in  Fig.  4  and  confirmed  by  the 
frequency  and  intrinsic  reaction  coordinate  (IRC)  calcula¬ 
tions  performed  at  the  CASSCF/AVTZ  level.  The  N-O  bond 
breaking  reaction,  which  leads  to  the  N(4S)+02(32g)  prod¬ 
ucts,  must  cross  over  a  high  barrier  of  1.15  eV  (TS1)  at  the 
single-point  MRCISD(<2)/AVQZ  level  of  theory.  In  contrast, 
a  small  barrier  of  0.03  eV  (TS2)  will  easily  result  in  the  0-0 
breaking  and  form  the  0('V)+N0(2H)  product  with  the  exo- 
thermicity  of  -AE(O+NO)_NOO=0.66  eV  at  the  same 
MRCISD(<2)/AVQZ  level.  Clearly,  the  NOO  species  pro¬ 
duced  via  the  CT  process  will  almost  exclusively  dissociate 
into  the  0(3P)+N0(2n)  products  due  to  the  extremely  small 
barrier  and  the  large  internal  energy  gained  from  the  CT 
reaction.  As  illustrated  in  Eq.  (4),  charge  transfer  will  be 
exothermic  for  neutral  reagents  X  with  an  IP<(IPad(NOO) 
-A£(0+no)_NOo)  =  10.44+0.66=  11.10  eV,  where  IPad(NOO) 
is  the  adiabatic  IP, 

NOO+  +  e-  — >  NOO,  AH  -  -  IPad(NOO)  =  -  10.44  eV, 


(4a) 

X-*X*  +  e~,  AH  =  IP x, 

(4b) 

N00->0(3P)  +  N0(2II), 

(4c) 

AH  =  A E(o+No)-noo  =  -  0-66  eV , 

NOO+  +  X  ->  X*  +  0(3P)  +  NO(2II) 

(4d) 

A  Htot  =  IP,  -  IPad(NOO)  +  A£(0+n0)  -NOO- 

Within  the  uncertainty,  the  energetics  agrees  well  with  the 
experimental  limit  of  1 1.2  eV.  In  essence,  the  overall  mecha¬ 
nism  is  dissociative  charge  transfer.  Typically,  dissociative 
charge  transfer  is  only  experimentally  observed  if  the  ionic 
product  dissociates,  which  is  easily  detected.  However,  the 
present  calculations  indicate  that  for  the  species  with  the 
highest  IPs  where  charge  transfer  occurred,  the  neutral  frag¬ 
ment  dissociation  must  also  have  occurred  and  might  also 
have  occurred  in  all  of  the  other  charge  transfer  reactions. 


C.  NOO++C2H6  potential  energy  surface 

While  detailed  calculations  for  all  of  the  reactions  in 
Table  I  are  not  practical,  we  have  performed  additional  cal¬ 
culations  on  the  hydride  abstraction  reaction  of  NOO+  with 
ethane  (C2H6).  The  optimized  reaction  intermediates,  transi¬ 
tion  states,  and  the  various  possible  products  are  shown  in 
Fig.  5  along  with  the  corresponding  potential  energy  profile. 
Due  to  the  size  of  the  system,  the  calculations  were  only 
carried  out  using  density  functional  theory  at  the  B3LYP/ 
aug-cc-pVTZ  level  of  theory.  The  product  HN02  structures 
in  Fig.  5  are  labeled  with  the  symmetry  of  the  species  in  the 
name  and  are  numbered  from  the  highest  energy  product 
starting  as  number  1  to  the  lowest  energy  product  of  a  given 
symmetry. 

Basically,  three  possible  reaction  channels  were  identi¬ 
fied  for  the  reaction  of  NOO+  with  C2H6.  The  two  lowest 
energy  channels  pass  through  transition  states  TS1  and  TS2, 
which  are  below  the  energy  of  the  reactants.  Both  TSs  lead  to 
C2H5  with  the  corresponding  HN02_CJ_  1  neutral  product, 
having  an  exothermicity  of -175.4  kJ  mol-1.  The  energies  of 
these  two  TSs  are  almost  the  same,  separated  by 
<2  kJmol-1.  Therefore,  they  might  contribute  similarly  to 
the  reactivity,  although  entropic  effects  would  also  play  a 
role.  The  fact  that  these  transition  states  have  lower  energies 
than  the  reactants  suggests  that  van  der  Waals  complexes 
must  exist.  These  complexes  are  labeled  as  INTI  and  INT2 
in  Fig.  5.  However,  van  der  Waals  complexes  were  not  lo¬ 
cated  at  the  current  level  of  computation.  This  must  be  re¬ 
lated  to  the  inability  of  the  density  functional  method  in  de¬ 
scribing  the  dispersion  interaction,  but  it  was  not  very 
important  for  understanding  the  overall  reactivity.  The  third 
possible  reaction  channel  would  proceed  through  the  higher 
energy  TS3  to  produce  an  isomer  of  HN02  as  the  neutral 
fragment  that  features  a  significantly  stretched  0-0  bond, 
given  as  HN02_Cf_4  in  Fig.  5.  Clearly,  the  major  difference 
between  the  third  channel  and  the  first  two  channels  is  that 
the  terminal  O  atom  of  NOO+  will  abstract  the  H  atom  from 
C2H6  along  the  highest  energy  pathway,  while  the  N  atom  of 
NOO+  will  abstract  the  H  atom  along  the  two  lower  energy 
pathways.  An  activation  energy  of  14.6  kJ  mol-1  would  be 
required  to  surmount  TS3,  which  makes  this  pathway  less 
important  at  low  energy  or  temperature.  Nevertheless,  reac- 
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FIG.  5.  ZPE  corrected  potential  energy  surface  for  the  reaction  of  NOO++C2H6.  Numbers  on  the  structures  refer  to  the  bond  length  (A)  and  bond  angles 
(degree).  The  vectors  on  the  TS  structures  correspond  to  the  reaction  coordinates  with  their  imaginary  frequencies  given  in  the  parentheses.  The  product  HN02 
structures  are  labeled  with  the  symmetry  of  the  species  in  the  name  and  are  numbered  from  the  highest  energy  product  starting  as  number  1  to  the  lowest 
energy  product  of  a  given  symmetry. 


lion  along  the  first  two  low  energy  pathways  could  proceed 
quite  easily,  which  explains  the  nearly  100%  efficiency  of 
this  reaction. 

V.  CONCLUSIONS 

The  NOO+  ion  was  found  to  be  highly  reactive  with 
most  neutral  molecules,  except  for  some  species  with  high 
ionization  potentials.  All  neutral  species  with  an  IP 
=£11.18  eV  reacted  with  NOO+,  in  part  by  charge  transfer. 
Detailed  calculations  showed  that  the  adiabatic  and  vertical 
IPs  of  NOO  were  about  1 0.4  and  1 1 .7  eV,  respectively.  The 
observed  experimental  limit  of  1 1 . 1 8  eV  agreed  well  with 
the  theoretical  charge  transfer  energetics  of  the  final  prod¬ 
ucts,  assuming  X++[0(2P)+N0(2n)]  formed  via  dissocia¬ 
tion  of  metastable  NOO  during  the  reaction.  The  neutral  re¬ 
agent  IP  must  thus  be  less  than  (IPd(NOO) 
-  A£(0+NO)_NOo)  =  ~  1  LI  eV  for  the  overall  charge  exchange 
reaction  to  be  exothermic.  This  case  would  be  dissociative 
charge  transfer  where  the  neutral  product  dissociated  instead 
of  the  ionic  product. 

Formation  of  NO+  product  ions  was  also  very  common. 
In  at  least  some  of  the  reactions,  this  product  ion  might  be 
due  to  collisions  of  NOO+  with  neutral  X  overcoming  a  dis¬ 
sociation  barrier,  especially  since  NOO+  is  separated  from 


NO++0  by  a  barrier  of  0.88  eV.6  Collision  induced  isomer¬ 
ization  to  the  lowest  energy  ONO+  isomer  was  also  found  in 
the  reactions  of  NOO+  with  CS2  and  S02.  Several  of  the 
reactions  had  numerous  product  channels.  The  CS2  reaction, 
in  particular,  resulted  in  six  different  products. 
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